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Abstract 

NMSPEC is a Fortran code that computes the sparticle and Higgs masses, as well 
as Higgs decay widths and couplings in the NMSSM, with soft SUSY breaking terms 
specified at Mqut- Exceptions are the soft singlet mass m? s and the singlet self coupling 
k, that are both determined in terms of the other parameters through the minimization 
equations of the Higgs potential. We present a first analysis of the NMSSM parameter 
space with universal SUSY breaking terms at Mqut ~ except for m s and A K - that 
passes present experimental constraints on sparticle and Higgs masses. We discuss in 
some detail a region in parameter space where a SM-like Higgs boson decays dominantly 
into two CP odd singlet-like Higgs states. 



PAC numbers: 12.60.Jv, 14.80. Cp, 14.80.Ly 



December 2006 
LPT Orsay 06-79 
LPTA Montpellier 06-61 

x Unite mixte de Recherche - CNRS - UMR 8627 
2 Unite mixte de Recherche - CNRS - UMR 5207 



1 



1 Introduction 



The Next to Minimal Supersymmetric Standard Model (NMSSM) [1,2] provides a very ele- 
gant solution to the \x problem of the MSSM via the introduction of a singlet superfield S. 
For the simplest possible scale invariant form of the superpotential, the scalar component 
of S acquires naturally a vacuum expectation value of the order of the SUSY breaking scale, 
giving rise to a value of // of order the electroweak scale. Hence the NMSSM is the sim- 
plest supersymmetric extension of the standard model in which the fundamental Lagrangian 
contains just SUSY breaking terms but no other parameters of the order of the electroweak 
scale. 

As in the MSSM, the phenomenology of the NMSSM depends on a certain number of 
parameters (mostly soft SUSY breaking parameters) that cannot be predicted from an un- 
derlying theory at present. It is then useful to have computer codes that compute physically 
relevant quantities as Higgs and sparticle masses, couplings, decay widths etc. as functions 
of the initial parameters in the Lagrangian. Such codes allow to investigate which regions 
in parameter space are in conflict with present constraints on physics beyond the standard 
model and, most importantly, which regions in parameter space can be tested in future expe- 
riments and/or astrophysical measurements. A review over corresponding publicly available 
computer codes can be found in [3]. 

Most of these codes are limited to the MSSM; the only presently available code for the 
NMSSM is NMHDECAY [4,5]. 

Many mechanisms for spontaneous SUSY breaking (such as supergravity with universal 
kinetic terms for all matter fields) predict a simple pattern for the soft SUSY breaking pa- 
rameters at a large scale as the GUT scale. Then, the soft SUSY breaking parameters are 
often universal, i.e. completely specified by universal gaugino masses, scalar masses squared 
and trilinear couplings. In the case of the MSSM, the corresponding scenario is denoted 
as the CMSSM (or mSUGRA since, in principle, other scenarios as GMSB or AMSB could 
constrain the soft SUSY breaking parameters). The advantage of such a scenario is, apart 
from its motivation through a simple underlying theory, obviously an enormous reduction 
of the number of unknown parameters. It would be desirable, if one could implement such 
simple boundary conditions also in the NMSSM, which we will denote as the CNMSSM. 

Numerical studies of the sparticle and Higgs spectrum in the CNMSSM have been per- 
formed in [2]. Since then, additional radiative corrections to the Higgs spectrum have been 
computed [6], and experimental constraints on the sparticle and notably the Higgs spectrum 
have become considerably stronger [7]. Also, the corresponding numerical codes have not 
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been made publicly available. 

NMSPEC is a code that allows to chose parameters for the NMSSM at the GUT scale. 
Exceptions are the singlet self coupling k and the soft SUSY breaking singlet mass m^; 
as described below, these two parameters are determined in terms of others through the 
minimization equations of the Higgs potential. Possible scenarios such as GMSB or AMSB 
are not (yet) implemented in NMSPEC. Also, like NMHDECAY, NMSPEC is at present limited 
to the simplest version of the NMSSM; neither terms linear in the singlet [8], additional U(l) 
gauge symmetries [9] or combinations thereof [10] are considered. (See [11] for a review of 
such more general versions of the NMSSM.) Subsequently we limit ourselves to a scale 
invariant superpotential of the form 

W = \SH u H d + ^S 3 + ... , (1.1) 

and the associated soft SUSY breaking trilinear couplings A\, A K and Higgs masses. 

In principle numerical codes, that are designed to compute the particle and sparticle 
spectrum depending on the soft SUSY breaking terms defined at the GUT scale, could proceed 
as follows: starting with numerical values for all unknown parameters at M GU t, the RG 
equations can be integrated numerically down to the weak scale. Then, the effective Higgs 
potential can be minimized numerically, and from the Higgs vevs at the minimum one obtains 
tan/3 and Mz (where Mz is used to determine a previously unspecified overall mass scale). 
This approach has been used for the NMSSM in [2]. 

It has several disadvantages, however: First, the numerical minimization of the Higgs 
potential is quite computer time consuming. Second, since the overall mass scale and tan f3 
are determined only at the end, one finds that m top , Higgs and sparticle masses often violate 
experimental bounds - it is not possible to input "large" soft terms (and the correct value 
of h t ) from the start. 

An advantage of this approach would be that it makes required (fine tuned) relations 
among the initial parameters immediately obvious. 

However, in view of its disadvantages, practically all present numerical codes (within the 
MSSM) proceed differently: both Mz and tan j3 are used as inputs, and the two minimization 
equations of the Higgs potential w.r.t. the two real Higgs vevs H u and are used to 
compute fi and the soft B parameter in terms of the others. (This leaves open the sign of //; 
subsequently the two possibilities have to be treated separately.) /i and B have only a small 
effect on the RG evolution of the other parameters, only via threshold effects from particles 
whose masses depend on \x and/or B. One can take care of this via a (typically rapidly 
converging) iterative procedure. (Now, however, possible fine tunings originating notably 
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from the requirement to chose B very precisely, are hidden.) 

At first sight, an application of this procedure to the NMSSM is not obvious: both \i 
and B are no longer independent parameters (although effective S dependent parameters 
/i e fj = XS, B cf[ = A\ + kS can still be defined), and one has to cope with three coupled 
minimization equations w.r.t. H u , Hd and S. 

A possible way out is the following: First, the tree level minimization equations w.r.t. 
H u and can be solved for fi e s and B e ^ , as in the MSSM, in terms of the other parameters 
(incl. Mz and tan/5). From /i e g- and B e g one can deduce (for A and A\ given) both S and 
k. Finally, from the minimization equation w.r.t. S, one can easily obtain the soft singlet 
mass m 2 in terms of all other parameters. 

This allows to chose as input parameters in the Higgs sector of the NMSSM the soft Higgs 
doublet masses squared, A\, A K , Mz, tan/5 and A. From these the three parameters /z c fr (or 
S), k and m 2 can be derived. The radiative corrections to the Higgs potential (that depend 
on all sparticle and Higgs masses) show a weak dependence on these parameters. They can 
be included in the minimization equations that become non-linear in the parameters to solve 
for, but one can solve them iteratively by a loop that converges rapidly. 

The derived parameters k and m 2 affect the RG evolution equations of some of the other 
parameters not only through threshold effects around M susy , but through the f3 functions 
(already at one loop). However, the numerical impact is relatively small such that an iter- 
ative procedure converges quite rapidly again. (These iterative procedures are described in 
section 2.) 

Clearly, once m 2 is an output rather than an input, it becomes difficult to find parameters 
such that m 2 at the GUT scale assumes the same value as, e.g., the Higgs doublet (or other 
scalar) soft masses squared. On the other hand it is easy to imagine that a mecanism 
for the generation of soft SUSY breaking terms treats the singlet differently from the other 
non-singlet matter multiplets. 

It happens frequently that the derived value of m 2 is negative even at the GUT scale. 
This does not imply that the singlet vev takes values of the order of the GUT scale, or that 
the scalar potential is unbounded from below: due to the term k 2 S a in the scalar potential, 
the singlet vev will always be of the order of m s j k (hence of the order of M susy ), and it is 
still necessary to integrate the RG equations also for the couplings and masses involving the 
singlet down to the SUSY scale and to compute further radiative corrections there. 

Once k is an output rather than an input, it is difficult to study the "Peccei-Quinn" limit 
k — > (with A finite) in the NMSSM [12]. However, this limit is quite unphysical at least 
in the simplest version of the NMSSM considered here: Once n tends to zero, the singlet 



4 



vev S tends to infinity, and hence n e g = XS becomes unacceptably large: large /j, e g implies 
a stable Higgs doublet potential at the origin (hence no electroweak symmetry breaking), 
unless at least one of the Higgs doublet SUSY breaking masses squared is of the same order 
and negative. Hence k — > would require M susy — > oo. 

The limit where k and A tend to zero simultaneously and the NMSSM turns into the 
MSSM plus a decoupled singlet sector poses no problem, however: it is easy to see that, for 
given and fixed /x e ff , B e g and Ax, A — > implies automatically k ~ A — »• 0. This limit is 
also stable under the RG flow from the GUT scale down to the electroweak scale. Although 
the vev S tends towards infinity in this limit, the masses of all components of the singlet 
superfield remain of the order of M susy (i.e. of the order of A K , Ax and fi c g). Hence, in order 
to study the MSSM limit of the NMSSM within NMSPEC, it is enough to chose a tiny value 
for A on input. (A = is not allowed, since expressions proportional to k/X would become 
ill defined. Also, A K still has to be chosen within a fi c g and Ax dependent window in order 
to avoid negative masses squared for the CP even and/or the CP odd singlet scalars.) 

Unless modified by the user, NMSPEC allows to chose the following input parameters: 
tan/3 and the sign of fi e g at the weak scale and A at the SUSY scale; at the GUT scale the 
free parameters are universal gaugino masses M1/2, universal scalar masses mg and universal 
trilinear couplings A . Exceptions are the following: 

a) the soft singlet mass at Mqut is an output, as described above; 

b) the trilinear coupling A K can be chosen differently from the other trilinear couplings. 
The reason for this is twofold: first, if an underlying mecanism for the generation of the 
soft SUSY breaking terms treats the singlet differently from the other matter fields (as it is 
already assumed for m 2 s ), this will also affect the coupling A K which involves singlets only. 
(We left Ax unified with the other trilinear couplings at Mgut; this could easily be changed 
by the user, however.) 

Also, interesting physics is associated with particular values of A K : in certain regions 
of the parameter space of the NMSSM, the lightest neutral CP even doublet-like Higgs 
boson (the one ressembling to the SM Higgs boson) can decay into two lighter neutral CP 
odd (sometimes also CP even) Higgs bosons [13,14]. This process allows for a doublet- 
like CP even Higgs boson h to escape LEP constraints even for nih below 114 GeV. (The 
corresponding LEP constraints have recently been updated in [7] , see the discussion in section 
3 below. The allowed window has become somewhat tighter, but still exists.) On the one 
hand, such scenarios can make the detection of a CP even Higgs boson at the LHC quite 
difficult [13, 14]. It has also been argued that the fine tuning, that is required in the MSSM 
in order to lift the Higgs mass above the LEP bound, is relieved in the corresponding region 
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of parameters in the NMSSM [14]. This region of parameters in the NMSSM requires, 
however, that A K assumes values within a certain range such that singlet-like Higgs bosons 
are sufficiently light. The possibility to chose A K differently from the other trilinear couplings 
in NMSPEC allows to study this range (now in terms of A K defined at M G ut), and the physics 
associated with it. 

In the next section 2 we describe in detail, how the code NMSPEC deals with the various 
input parameters in the NMSSM. In section 3 we discuss new features of some of the subrou- 
tines of NMHDECAY v3; most of the subroutines of NMHDECAY v3 are shared by NMSPEC. 
In section 4 we present first results obtained with NMSPEC. The allowed parameter space 
(passing all phenomenological tests) divides into three regions: a MSSM-like region with A 
relatively small, and two NMSSM specific regions with A large either for small tan/5 and 
above 114 GeV, or for medium to large tan f3 and rrih below 114 GeV, but A K within a certain 
window that renders the unconventional Higgs decay modes discussed above possible. The 
width of the corresponding window (and the fine tuning associated with it) is discussed in 
some detail. 

2 The code NMSPEC 

As discussed in the introduction, the code NMSPEC proceeds as follows: Input parameters 
are the SM gauge couplings, Fermion masses of the 3rd generation, tan/3, the sign of /x e g, 
and A as well as the soft SUSY breaking masses and trilinear couplings (with the exception 
of m 2 s ) at the GUT scale. 

The code NMSPEC proceeds as follows: First, a guess for k, m 2 , n e g and M susy is made, 
and a preliminary calculation of Mgut and the gauge and Yukawa couplings at the GUT 
scale is performed through an integration of the (two loop) RG equations from the SUSY 
scale up to Mgut using the subroutine RGES. (Mgut is defined as the scale where ol\ and 
a 2 unify.) 

Then, an "external" loop starts, that involves both the GUT scale and M susy : the (two 
loop) RG equations for the soft terms, together with the gauge and Yukawa couplings, are 
integrated from Mgut down to M susy in the subroutine RGESINV. 

There an "internal" loop starts, that involves physics between M susy and the weak scale: 
In the subroutine RUNPAR, first the mass M\ of a potentially heavy Higgs doublet is 
estimated (which is required for the threshold effects below). The SUSY scale (denoted by 
Q2) is defined in terms of the squark masses of the first generations, and a scale QSTSB 
is defined in terms of the left and right stop masses squared. Then, the SM gauge and 
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Yukawa couplings are calculated at the scale QSTSB (together with the Higgs wave function 
renormalization constants Zj). Threshold effects at Higgs and sparticle masses are taken care 
of. Finally, the NMSSM specific Yukawa and trilinear couplings are computed at QSTSB. 

In the next subroutine MSFERM, the (/z e ff-dependend) 3rd generation squark and slepton 
pole masses and mixing angles are computed at the scale QSTSB. They are needed below 
for the radiative corrections to the minimization equations of the Higgs potential. 

In the following subroutine LOWMUK, first the RG equations for the soft Higgs doublet 
masses squared are integrated from Q2 to QSTSB. Then, /j, e g, k and m 2 are determined 
from the minimization equations of the Higgs potential including radiative corrections, that 
depend on previous estimates for /x e g and k. A variable CHECK verifies the relative change of 
/i c fj and k with respect to its previous values. Finally, an improved value for M\ is computed 
(in case the internal loop stops here; otherwise its previous estimate is used again). 

Unless CHECK < 10~ 6 , the internal loop returns now to the subroutine RUNPAR. There, 
the improved values for // e g and k are used for improved threshold corrections. In the next 
call of LOWMUK, the radiative corrections to the minimization equations also include the 
improved values of /i e s and k. 

If CHECK < 10~ 6 , fi c s and k, have been determined to a sufficient accuracy. We have 
never observed a convergence problem related to the determination of /x e ff and k; 5 inner 
loops are usually more than sufficient. 

Then, the "external" loop continues: next, the RG equations for the SM gauge and 
Yukawa couplings are integrated again from the weak scale up to Mqut i n the subroutine 
RGES, using improved threshold effects at sparticle and Higgs masses as well as the updated 
value for k. The value for Mqut is improved simultaneously. 

Next, the RG equations for the soft terms are integrated from M susy up to Mqut in the 
subroutine RGESUNI, using the new SM gauge and Yukawa couplings and the new values 
for Mqut and m 2 . A variable GUTEST verifies, whether all soft terms coincide with the 
required boundary conditions at Mqut- 

Unless the relative deviation determined by GUTEST is less than 10~ 4 (for all soft SUSY 
breaking parameters), the external loop jumps back to a call of RGESINV. There, the new 
values for Mgut and the SM gauge and Yukawa couplings are used in order to integrate 
again the RG equations for the soft terms (with the required boundary conditions at Mqut) 
down from the GUT scale. 

If GUTEST is less than 10~ 4 , the required values of the soft terms at Mqut are considered 
as sufficiently precise, and the external loop is left. The subsequent subroutines of N MS PEC 
are shared with NMHDECAY v3: Higgs, gluino, chargino and neutralino masses are computed, 
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and all Higgs branching ratios into SM particles and sparticles are determined. Then, present 
experimental constraints on Higgs and sparticles masses and couplings are verified. Warnings 
indicate at the end, whether such constraints are violated. A rough (1 loop) computation 
of the b — > s 7 branching ratio is performed (for information only, without generating an 
error message), and it is verified whether the choice of parameters corresponds to a global 
minimum of the NMSSM Higgs potential. 

Finally, if desired (if the flag OMGFLAG is on), the dark matter relic density is com- 
puted using a NMSSM version of the micrOMEGAs_2.0 subroutines [15]. These subroutines 
are called internally in NMHDECAY v3 and NMSPEC; an interface file, as it was the case for 
NMHDECAY v2 [16], is not needed anymore. The NMSSM version of micrOMEGAs_2.0 is in- 
cluded in the downloadable package NMSSMTools_1.0, and compiled together with NMHDE- 
CAY and NMSPEC, as described in the appendix. 

The different options (a single point in parameter space, scans over regions in parameter 
space using grids or random numbers) are also described in the appendix. 

3 New features in NMHDECAY v3 and NMSPEC 

The code NMHDECAY v3 performs the same tasks as the codes NMHDECAY vl and v2, that 
are described in detail in [4,5]. Now, many subroutines are shared by NMSPEC, and the 
codes NMHDECAY v3 and NMSPEC can be downloaded simultaneously (see the appendix). 
Apart from the architecture, some properties of some subroutines have been improved with 
respect to NMHDECAY v2; these changes are described below. 

3.1 Experimental constraints 

The (negative) results of the four LEP collaborations ALEPH, DELPHI, L3 and OPAL on 
Higgs Boson searches have been updated and combined by the LEP Higgs Working group 
in [7]. As compared to the corresponding constraints already implemented in NMHDECAY 
v2, the updates in [7] include improved constraints on Higgs masses vs. couplings for the 
following processes: 

• Z — > ZH 2 — > ZH 1 H 1 where H 1 stands for a CP even or a CP odd neutral Higgs boson. 
The considered decay channels of Hi are Hi — » 2b and Hi — > 2r, and combinations 
thereof. 

• Z — > H1H2 where both Hi and Hi can denote a CP even or a CP odd neutral Higgs 
boson. Decays of both Hi and Hi into 2b and 2r are studied. 
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• Z — > H1H2 — > H\H\H\. Only simultaneous decays of all three ifi states either into 66 
or into 6r have been considered. 

The corresponding data files are included in the directory EXPCON. They restrict the 
allowed parameter space of the NMSSM more strongly than within the previous versions vl 
and v2 of NMH DECAY. In particular, the constraints on the decay of a SM-like Higgs boson 
h to aa, where a is a CP odd (singlet-like) scalar, are much stronger. Roughly speaking, the 
only remaining possibilities (unless the Z-Z-h coupling is heavily suppressed) are: 

• 86 GeV < m h < 106 GeV and m a < 11 GeV (where the bb decay of a is impossible, 
and no constraints from a — > tt are available), 

• 106 GeV < m h < 114 GeV with m a going up to m h /2. 

(It must be noted, however, that LEP constraints are implemented only for individual pro- 
cesses. Combinations of different processes can rule out points in parameter space, that seem 
to have passed the individual constraints.) 

3.2 Input parameters 

The implicit definitions of the input parameters A, k, /z eff , A x and A K in NMHDECAY v3 differ 
slightly from the ones in NMHDECAY v2: now these are all defined at a common SUSY scale 
Q2 (essentially the average squark masses of the first two generations, unless Q2 is defined 
by the user), in the DR scheme as all other soft SUSY breaking parameters. The numerical 
effect concerns essentially A\. in order to reproduce data points obtained with the previous 
versions of NMHDECAY, A x has to be shifted by an A top dependend amount (that can be 
estimated from the RG equation for A\). 

3.3 Higgs mass 

The calculation of the radiative corrections to the Higgs masses has been reorganised and 
improved: Now all parameters used by the corresponding subroutine MHIGGS (gauge and 
Yukawa couplings, tan /3, trilinear couplings, quark and squark masses, fi c g and Higgs vevs) 
are defined at a scale QSTSB, that corresponds to an average of the right and left stop 
quark masses. (QSTSB is computed separately from the scale Q2 above, and can be lower.) 
Corresponding formulas for the radiative corrections have essentially become simpler, except 
that contributions from gauge boson and light Higgs boson loops have now to be added 
explicitly. (In the LLA, which is used for the elektroweak corrections, contributions from 
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particles with masses ~ Mz to masses of particles of the same order do not involve large 
logarithms. Now, the results on heavy Higgs masses are more reliable, without negative 
effects on the precision of the light Higgs Boson mass.) 

Inside the subroutine MHIGGS, dominant gaugino loop contributions are now expressed 
in the form of a shift of A\, which also improves the precision of heavy Higgs spectrum. 

New radiative corrections, that have not been considered before, are two loop contribu- 
tions to the lightest Higgs mass involving two powers of large logarithms and electroweak 
gauge couplings gi/2 of the form ~ h\ g 2 j 2 log 2 , where log denotes a potentially large log- 
arithm. The dominant effect of these contributions depends on the chargino/neutralino 
spectrum, and is described in [17]. 

4 The allowed parameter space of the CNMSSM 

In the CMSSM, the present experimental constraints on sparticle and Higgs masses can be 
satisfied by chosing sufficiently large soft SUSY breaking terms at the GUT scale, and tan (3 
not too small. As stated above, the required fine tuning is hidden through the implicit 
determination of /i and B. 

In the CNMSSM (without unification of the soft singlet mass, see above) it is often nec- 
essary, in addition, to chose a non-universal value for A K in order obtain positive masses 
squared for both CP even and CP odd singlet scalars, even when they are essentially decou- 
pled. Otherwise, for A < 0.15, the allowed parameter space of the CNMSSM corresponds 
essentially to the one of the CMSSM (except for the possibility to have below 114 GeV 
and unconventional Higgs decay modes, see below). 

Larger values for A are not always allowed within the CNMSSM, however; possible prob- 
lems arise typically in the CP even Higgs sector: on the one hand, non- vanishing A increases 
the mass of the lightest Higgs doublet, at tree level, by an amount 

Ami = 2\ 2 M 2 z sm 2 2[3/(gl + g 2 2 ) . (4.1) 

This increase is relevant, however, only for small tan/3. On the other hand non- vanishing 
A induces mixings between the doublet and the singlet states. Since the singlet state is 
typically quite heavy, this mixing reduces the mass of the lightest eigenstate of the 3x3 CP 
even mass matrix. This lightest eigenstate will then often violate bounds from LEP [7], in 
spite of its singlet component. 

In order to maximize the mass of the lightest CP even mass eigenstate, this singlet- 
doublet mixing has to vanish (as described in [17]), which implies a relation between /i e fj, 
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Figure 1: Allowed values for A vs. tan/5 for M 1/2 = m = 500 GeV, A = -800 GeV, 
A K = -1500 GeV. 

A\, A, tan (3 and k. This relation is generally not satisfied within the CNMSSM; then - at 
least for larger values of tan (3 - this mixing effect disallows often large values of A. 

For given soft SUSY breaking terms at the GUT scale, one then obtains typically allowed 
"islands" in the plane A vs. tan (3. An example is given in the Figure 1 below, which shows 
such an allowed island for M 1/2 = m = 500 GeV, A = -800 GeV, A K = -1500 GeV and 
/Ucfr > (with mtop =175 GeV). 

The lightest (doublet-like) Higgs mass is ~ 114 GeV near both ends of this region, and 
up to ~ 117.5 GeV near its center. All three Yukawa couplings A, k and h t are relatively 
large, and the region is bounded - apart from LEP constraints on - by the condition on 
the absence of a Landau singularity in the Yukawa sector below Mgut- 

Larger values of A are also possible for larger values of tan /3. Then, however, the light- 
est (still doublet-like) Higgs mass falls below 114 GeV, and A K has to be chosen within a 
relatively narrow window in order to generate a light (singlet-like) CP odd scalar into which 
h can decay [13, 14]. Such a line of allowed values for A K vs. A is shown in Figure 2, for 
tan/5 = 5, Mi/a = m = 500 GeV and A = -800 GeV. 
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Figure 2: Allowed values for A K vs. A for tan/3 = 5, My 2 — m — 500 GeV and A = 
-800 GeV. 

For A < 0.16, this line joins the CMSSM-like region where rrih > 114 GeV. Along the 
line, rrih decreases to rrih ~ 86 GeV for A ~ .4. Larger values of A, i.e. lower values of rrih, 
are excluded by LEP. Across the line, the mass of the singlet-like CP odd scalar goes up to 
rrih/2 for A < .27 (where rrih > 106 GeV), but m a < 11 GeV for larger values of A (where 
m h ^ 106 GeV); a larger value for m a would imply that a decays mainly into bb, which is 
excluded by LEP. The width of the line, expressed in terms of A K , is ~ 1.7 GeV for A < .27, 
but shrinks down to ~ 100 MeV for larger values of A. This allows to estimate the fine 
tuning in A K required for this region of the parameter space of the CNMSSM. 

In principle, the width of the line (and hence the required relative fine tuning) can depend 
on all other soft SUSY breaking terms. In order to see this, we have chosen smaller values for 
Mi/a = 200 GeV, m = 100 GeV, fixed A = 0.2 (still at tan/3 = 5), and plotted the allowed 
values for A K vs. A in Figure 3. 

In Figure 3, rrih ~ 90 GeV near Aq ~ —200 GeV and increases up to ~ 112 GeV for 
A - -700 GeV. (For values of A < 700 GeV one finds m Ti < 100 GeV, which is excluded 
by stop searches at LEP). The mass of the singlet-like CP odd scalar is again < 11 GeV for 
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Figure 3: Allowed values for A K vs. A for tan/? = 5, My 2 = 200 GeV, m = 100 GeV and 
A = 0.2. 

A > - 440 GeV, but varies from a few GeV up to m h /2 for A K < - 440 GeV. The width 
of the line, expressed in terms of A K , is now ~ 2.5 GeV for A ~ —700 GeV, i.e. the relative 
fine tuning in A K is somewhat less than before. According to our preliminary analysis, this 
width cannot be increased by chosing different values of tan (3. 

On the other hand allowed regions in parameter space with < 1 14 GeV remain present 
also for smaller values of A (at least down to A ~ 0.1). The present analysis of the allowed 
parameter space of the CNMSSM is far from complete, but NMSPEC will allow for more 
detailed studies in this direction in the future. In any case it is important to note that the 
region in the parameter space of the NMSSM where the lightest SM-like Higgs boson decays 
dominantly into two CP odd scalars remains present in the CNMSSM; we recall that it could 
then be quite challenging to detect just one Higgs boson at the LHC [13]. 
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Appendix: How to use NMSSMTools 



The package NMSSMTools_1.0 contains NMSPEC, the new version of NMHDECAY (v3), a 
NMSSM version of micrOMEGAs_2.0 and an updated version of the LEP constraints. It can 
be downloaded from the web page: http://www.th.u-psud.fr/NMHDECAY/nmssmtools.html. 

Once this compressed tar file is downloaded, one has to type 

tar -zxvf NMSSMTools. 1 . tgz 

which will create the directory NMSSMTools_1.0. Inside this directory one finds: 

• the directory sources, that contains the subroutines common to NMHDECAY and 
NMSPEC, the directory micromegas_2.0.3 with the subroutines for the computation 
of the dark matter relic density, and the directory EXPCON with data files used for 
the experimental constraints; 

• the main programs for NMHDECAY and NMSPEC: nmhdecay.f, nmhdecay_rand.f, 
nmspec_grid.f, nmspec.f, nmspec_rand.f and nmspec_grid.f; 

• six sample input files: inp.dat, randinp.dat, gridinp.dat, inpsp.dat, randinpsp.dat and 
gridinpsp.dat. All input (and output) files are in SLHA format [18], including some 
NMSSM specific switches that are described in [4,5,19]. We have added privately 
defined switches for the boundary values of the parameters to be scanned over, the 
number of points to be scanned over, and an integer ISEED that serves as input for 
the random number generator. A sample input file randinpsp.dat is given in Table 1; 

• corresponding output files: spectr.dat, decay.dat, omega. out, randerr.dat, randout.dat, 
griderr.dat, gridout.dat, spectrsp.dat, decaysp.dat, omegasp.dat, outsp.dat, randerrsp.dat, 
randoutsp.dat, griderrsp.dat, gridoutsp.dat; 

• A master Makefile. 

In order to compile the subroutines in the directories sources and micromegas_2.0.3, one has 
to type first: 

make init 

Then to compile the main programs, type: 
make 
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This will create 6 executable files: nmhdecay, nmhdecay_rand, nmhdecay_grid, nmspec, 
nmspec_rand and nmspec_grid. Now the user has the following options: 

• To study a single point in parameter space, with input parameters at the SUSY 
scale: Edit the input file inp.dat. Type nmhdecay. The output spectrum is in the 
file spectr.dat, the branching ratios of the 6 Higgs states in the NMSSM are in the 
file decay.dat. The dark matter relic density (if OMGFLAG=l in inp.dat) is in the file 
omega.dat. (A first call of the computation of the dark matter relic density provokes 
the compilation of additional subroutines of micrOMEGAs_2.0.) 

• To scan over randomly chosen points within specified boundary values, with input pa- 
rameters at the SUSY scale: Edit the input file randinp.dat. Type nmhdecay_rand. The 
number of points that have passed all tests and their ranges are in the file randerr.dat. 
The output is in the file randout.dat. In order to specify the output, you need to edit 
the file nmhdecay _rand.f which contains the subroutine OUTPUT. At the beginning of 
this subroutine, you find the different elements of the array RES, whose elements are 
printed into the output file randout.dat. You can declare the elements of the array RES 
according to your needs; note that the input parameters are elements of the array PAR 
whose meaning is given at the beginning of the MAIN program in nmhdecay _rand.f. 
After saving a modified version of nmhdecay_rand.f it is necessary to type make again. 

• To scan over a lattice of points within specified boundary values, with input parameters 
at the SUSY scale: Edit the input file gridinp.dat. Type nmhdecay ^grid. The number 
of points that have passed all tests and their ranges are in the file griderr.dat. The 
output is in the file gridout.dat. In order to specify the output, you need to edit the 
subroutine OUTPUT in the file nmhdecay _grid.f, and proceed as described in the case 
of a random scan above. 

• To study a single point in parameter space, with input parameters at the GUT scale: 
Edit the input file inpsp.dat. Type nmspec. The output spectrum is in the file spec- 
trsp.dat, the branching ratios of the 6 Higgs states in the NMSSM are in the file 
decaysp.dat. The dark matter relic density (if OMGFLAG=l in inp.dat) is in the file 
omegasp.dat. The SUSY scale parameters are written in the file outsp.dat in the same 
format as the file inp.dat, so that outsp.dat can be used as an input file for nmhdecay. 

• To scan over randomly chosen points within specified boundary values, with input 
parameters at the GUT scale: Edit the input file randinpsp.dat (cf table 1). Type 
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nmspec_rand. The number of points that have passed all tests and their ranges are 
in the file randerrsp.dat. The output is in the file randoutsp.dat. In order to specify 
the output, you need to edit the subroutine OUTPUT in the file nmspec_rand.f, and 

proceed as described above. 

• To scan over a lattice of points within specified boundary values, with input parameters 
at the GUT scale: Edit the input file gridinpsp.dat. Type nmspec_grid. The number 
of points that have passed all tests and their ranges are in the file griderrsp.dat. The 
output is in the file gridoutsp.dat. In order to specify the output, you need to edit the 
subroutine OUTPUT in the file nmspec_grid.f, and proceed as described above. 

To delete all the already compiled codes, the user should type: 
make clean 
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# INPUT FILE FOR NMSPEC VERSION 1 

# BASED ON SUSY LES HOUCHES ACCORD II 



RT OCK 


MODSFT 

1 II ll/lll ll ■ 






Q 
o 


J. 


it 

it 




RT ncv 


DllXlM.rU ID 






1 


127.92D0 


# 


ALPHA_EM--1(MZ) 


2 


1.16639D-5 


# 


GF 


■5 

o 


011 79nn 




AT PHA SfM/") 


4 


Q1 1 R7no 

c7 J. . lOl UU 


# 


MZ 


c. 


a 91 ado 


# 


MB (MB), RUNNING B QUARK 


6 


175. DO 


# 


TOP QUARK POLE MASS 


7 


1 . 777D0 


# 


MTAU 


BLOCK 


MINPAR 









0. 


# 


REN . SCALE 


4 


1. 


# 


SIGMU 


17 


100. 


# 


M0MIN 


18 


100. 


# 


M0MAX 


27 


200. 


# 


M12MIN 


28 


200. 


# 


M12MAX 


37 


5.0 


# 


TBMIN 


38 


5.0 


# 


TBMAX 


57 


-700. 


# 


AOMIN 


58 


-550. 


# 


A0MAX 


BLOCK 


EXTPAR 






617 


.2D0 


# 


LMIN 


618 


.2D0 


# 


LMAX 


647 


-80. 


# 


A_KMIN 


648 


-55. 


# 


A_KMAX 


BLOCK 


STEPS 









1000 


# NT0T 


1 


-1 


# 


I SEED 



Table 1: The randinpsp.dat file for sample parameter scan. 
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